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Electric current fluctuations are one type of unavoidable machine imperfections and induce magnetic-field 
perturbation as a source of instabilities in accelerators. This paper presents the measurement-based modeling 
methodology of fluctuating electric current arising from the power system of Fermilab's Booster synchrotron to 
discuss the ramifications of the presence of ripple current and space-charge defocusing effects. We also present 
the method of generating stochastic noise and the measurement and analysis methods of ripple current and 
offending electromagnetic interferences residing in the Booster power system. This stochastic noise model, ac- 
companied by a suite of beam diagnostic calculations, manifests that the fluctuating power-supply current, when 
coupled to space charge and impinging upon a beam, can substantially enhance beam degradation phenomena — 
such as emittance growth and halo formation — during the Booster injection period. With idealized and uniform 
charge-density distribution, the fractional growth of rms emittances due to ripple current under space charge turn 
out to be about 8 ~ 9 % in both transverse planes over the injection period of 2.2 ms prior to beam acceleration. 
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I. INTRODUCTION AND MOTIVATION 

As is common in other natural systems, subtle fluctuations 
are ubiquitous and inevitable in particle accelerator systems. 
An ensemble of charged particles is defined as a system, and 
all the beamline components (magnets, power supplies, RF 
cavities, beam position monitors, etc.) for accelerating, guid- 
ing, and diagnosing particle beams as environment, or sur- 
roundings. The system of a charged-particle beam perceives 
the environment of beamline components as a source of noise 
as illustrated by FIG. \T\ External noise is intrinsic to parti- 
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FIG. 1: System and surroundings 

cle accelerators of all types due to unavoidable machine im- 
perfections; e.g., ripple current from power supplies, ground 
vibration motion, etc. After being motivated by earlier find- 
ings from an idealized and simplified theoretical model 0], 
we speculated that the adverse influence of power- supply cur- 
rent fluctuations possibly account for beam loss phenomena 
observed during the injection process of the Booster. Hence, 
an independent and novel approach has been made to build a 



realistic stochastic noise model, based upon a series of power- 
supply noise 1 measurements, to investigate and explore the 
impact of current fluctuations on charged-particle beams of 
the Booster at injection energy of 400 MeV. 



II. MODELING METHODOLOGY 

For investigations of the impact of power- supply current 
fluctuations on charged-particle beams in an accelerator lat- 
tice structure, we began with building a preliminary noise 
model. As a first step, we designed and added a new physics 
module for generating stochastic noise to the existing ORBIT- 
FNAL package H. The new noise module is capable of gen- 
erating a wide spectrum of stochastic noise employing the 
Ornstein-Uhlenebck stochastic process [ 3] that is governed by 
a Langevin-like stochastic differential equation Jj] . 

Prior to detailed experimental measurements, we corrob- 
orated with the preliminary noise model using a linear lat- 
tice that non- white, or colored noise could possibly enhance 
beam degradation process of our interest. These preliminary 
findings are consistent with earlier findings from the theoreti- 
cal model of collective space-charge modes coupled with dy- 
namic noise 01]. As a next step, in addition to adding the new 
noise module to the ORBIT-FNAL, the existing TeaPot mod- 
ule was upgraded to establish Fermilab's Booster ring using a 
TEAPOT-style | 5] Booster lattice, and the Diagnostic module 
with new parallelized beam diagnostic calculations: actions, 
halo amplitudes, etc. 

As confirmed with the preliminary model, we proceeded 
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In the present context, we will use the terms (power-supply) noise and 
current fluctuations and ripple current interchangeably. 
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FIG. 2: Multi-stage approach to modeling the impact of power- supply noise on a beam with realism 




FIG. 3: Mapping from an experimental signature space to a 
stochastic parameter space 



A. External Noise 

External noise is originated from a source outside of the 
system; that is, beamline components. Since the effects of ex- 
ternal noise can be described by a stochastic differential equa- 
tion (SDE) 2 , we chose Langevin Equation (LE) as SDE for 
indeterministic current fluctuations arising from each GMPS 
unit 3 . It should be noticed that power- supply ripple currents 
are considered indeterministic, or random, or aperiodic, in the 
sense that it never exactly repeats itself. 



to devise methods for direct measurements of common-mode 
and differential-mode voltages, and ripples in the electric cur- 
rent. Repeated measurements and Fourier analysis confirm 
that a substantial amount of noise, which can be transmit- 
ted to the magnet system, is indeed present in the power sys- 
tem. Moreover, we performed equivalent-circuit simulations 
to investigate any offending resonances that can float around 
the magnet system. Based upon the measurement data and 
the results of Fourier analysis, stochastic parameterization of 
Booster ripple current is performed by means of matching 
power spectral densities between measured ripple currents and 
modeled Ornstein-Uhlenbeck (O-U) noise | 3]. While translat- 
ing modeled O-U noise to induced magnetic-field fluctuations, 
we tracked macroparticles in the presence of 3-D space-charge 
effects. 

The following FIG. [2] and FIG. [3] illustrate the multi-stage 
approach to the stochastic noise modeling. For the purpose 
of the stochastic parameterization of ripple current, mapping 
from an experimental signature space to a stochastic param- 
eter space was done; the preliminary model was faithfully 
tuned up with the power- supply noise measurements. As a 
consequence, we managed to match FFT power spectral den- 
sities between physical noise and modeled O-U noise. 



III. CLASSIFICATIONS OF NOISE 



B. Electromagnetic Interferences 

Electromagnetic-Interference (EMI) noise 10, \^ results 
from rapid changes in voltage and current in a power sup- 
ply. Transmissions of EMI noise are characterized as ei- 
ther radiative, or conductive. Conductive EMI noise, such 
as differential-mode (DM) and common-mode (CM) noise, 
is usually several orders of magnitude higher than the radia- 
tive EMI, and can be more harmful to the system. Given 
impedance (Z(co)) as a function of frequency CO, fluctuations 
in common-mode voltage (Vcm) create common-mode current 
(Icm), in addition to the inherent ripple current arising from 
sudden potential changes in the power-supply system. The 
EMI problem is thereby worsened and could result in larger 
current fluctuations, or common-mode current, or severe sys- 
tem damage. 



IV. STOCHASTIC PROCESS 

Of several different stochastic processes 4 , such as Pois- 
son process, Wiener process, etc., we chose the Ornstein- 
Uhlenbeck process to represent electrical current fluctuations 
as classified in the preceding subsection. 



In general, noise can be categorized into two types: external 
noise and internal noise. In our stochastic noise model, ripple 
current arising from Gradient Magnet Power Supply (GMPS) 
units of the Booster are regarded to be fluctuating external 
influences acting on the Booster beam (see FIG.Q]). 



2 A stochastic differential equation is a differential equation (DE) in which 
the coefficients are random functions of the independent variable, which is 
usually time. 

3 The Booster power supply will be frequently referred to as GMPS. 

4 In this paper, we use terms stochastic process, stochastic function, random 
process, random function and process interchangeably. 
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A. Langevin Equation 

In 1908, after the formulation of the Brownian movement 
by Einstein and Smoluchowski tlOU, P. Langevin introduced 
the concept of the equation of motion of a stochastic vari- 
able (i.e., the position coordinate of a Brownian particle) fTTh . 
Langevin Equation (LE) is considered to be the first example 
of a Stochastic Differential Equation (SDE) 5 . Langevin wrote 
the equation of motion for a Brownian particle according to 
Newton's second law under the assumption that a Brownian 
particle is subject to two forces: damping force (fd) and fluc- 
tuating force (ff) 



F(0 = Td{t)+r f {t) 



(1) 



d\(t) 
dt 2 



dt 



+ 



systematic force: dissipation 



stochastic force: fluctuation 

(2) 

Here, F (t ) , m , * , and £ represent the total force, particle mass, 
displacement, and the friction coefficient, respectively. The 
first term (— £ ?t(t) 6 ) on the right-hand side of Eqn. (|2]) rep- 
resents the viscous drag as a function of time, or dynamic 
friction. The second term ff(t) represents fluctuations which 
could be from white noise 7 or non- white noise. The form of 
Eqn. (0 can be transformed to Eqn. ([3]) of first order. For 
modeling physical noise of the Booster power system, we em- 
ployed non-white noise, or off-white noise in our investiga- 
tions. 



v(t)+av(t) =£(t), 



(3) 



where L (t) is a stochastic driving force, and a represents tjm. 
The following assumptions are made about the fluctuation part 

L(t): 

(1) L (t) is a function of time only, and independent of ^. 

(2) zero-mean distribution; 

(L(t))=0 (4) 

(3) The variation rate of L (t) is much faster than the velocity 
of a Brownian particle, v(t). Hence, the autocorrelation 
function Cl (t, t') is of Dirac-8 function; i.e., 



{L(t)L(t'))=Ah(t-t') 



(5) 



The expressions above define the statistical properties of 

L (t). There is a great advantage in using LE instead of us- 
ing Fokker-Planck Equation (FPE) of mathematical complex. 
The Langevin 's method is much easier to understand than the 



FPE since it is based upon the time evolution of a stochastic 
variable, whereas the FPE applies to the time evolution of the 
probability distribution. As such, LE allows us to dispense 
with the calculation of the diffusion coefficient, thus reducing 
associated mathematical complications. As a consequence, 
we built an effective but far more simplified model of stochas- 
tic noise. 



B. Ornstein-Uhlenbeck Process 



As in Eqn. 0]), LE for Brownian motion is given as, 



4(0+g«(0 = ^(0 



(6) 



As explained in the preceding subsection, LE is associated 
with 8-correlated Gaussian stochastic forces of statistical 
properties. This stochastic process 5(0 is referred to as an 
O-U process. Note that the noise strength A does not depen- 
dent upon the variables 5- Since the O-U stochastic process 
is inherently to represent the velocity of a Brownian particle, 
it is the appropriate choice of stochastic process for modeling 
electric noise, or current fluctuations (3|] ; i.e., a time-derivative 
of electric charge {dQ/ dt). Both position (x(t)) and velocity 
(v(0) describe Langevin's Brownian motion. However, by 
utilizing the O-U process of the velocity of a Brownian parti- 
cle, LE can be reduced to a I st -order linear stochastic differ- 
ential equation that is derived from Newton's 2 nd law. As a 
result, LE as a I st -order SDE is straightforward to find solu- 
tions. 

The O-U process is associated with an exponentially- 
decreasing autocorrelation function C$(t, t') lfl2ll and a finite 
autocorrelation time x ac JH] 8 . 



Cg(f, = (5(05(0 > =^exp(-(O flC |f - 0), 



where l n 



(7) 



in which 5(0' and A are a stochastic function, an auto- 
correlation frequency, and a constant noise strength, respec- 
tively. O-U processes are associated with the following prob- 
ability density function: 



^i(5, = 



i 



y/nJl/(i} a 



• exp 



\ 2A/2a) ac J 



(8) 



According to the Doob's theorem [13], the O-U process is the 
only stochastic process with all of the following properties: 
(1) stationary process, (2) Gaussian process, (3) Markovian 
process. In particular, if a process is invariant to translations 
in time (e.g. a shift in time (a)) then the process is called a 
stationary process 0, [3. For a stationary process, we can 



5 A stochastic differential equation is a differential equation with a stochas- 
tic (random) term. Therefore, its solution is also a random function. 

6 The notations, i and 'i , denote d^/dt and d 2 * / dt 1 , respectively. 

7 White noise is noise with a flat frequency spectrum. 



The autocorrelation function c(t, t') determines the property of a stochas- 
tic process of interest. 
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make the following simplifications: 

<5(*i + a)%(t 2 + a) ••• %[t n + a) ) = (£('1)^2) ••• $(t n ) ) 

(9) 

where b,(t) is the stochastic function and ( ... ) is the statisti- 
cal average. The form of Eqn T d9l) implies the fallowings • 

(1) ^lnce^the underlying mechamsms causing the fluctuations 

do not change with time, the stochastic properties of a sta- 
tionary process are conserved. 

(2) The important parameter in the O-U process is relative 
time and not the absolute time. 

(Uh-aMh)) = a(hMh + «)) (10) 

Therefore, the autocorrelation function C^(t, t') for a sta- 
tionary process is a function of \t — 1'\ only. 

(3) The ensemble average and the time average are the same, 
which leads to the ergodic property. 

C. Markov Process 

Markov processes involve the use of conditional probabil- 
ity. 

w 2 £u & t) = ^1(^1)^2(^1 1 £2, (ii) 

The Markov process is therefore defined as follows 0] : 

$n( £l t\, ^2t2, • • • , Z>n-1 t n -\ \ t>n t n ) = Z>n-1 t n -\ \ £>nt n ) 

(12) 

The form of Eqn. (IT2l) implies that all the ^ n for n > 2 can be 
derived, when only ^2 is known. In other words, 

Only the present condition determines the future condition. 

In order to avoid unnecessary mathematical complexity in 
building a stochastic noise model, we exploited the Marko- 
vian property. This is another reason why, of all the stochastic 
processes, the O-U process is the most convenient choice for 
modeling the Booster ripple current. 

V. NON-WHITE NOISE GENERATION 



tion function Cx\ : 

c„(f, t') = (Ti(f)T!(0) = ^-e X p(-co ac |f-r'|) (14) 

Non- white noise T| is governed by LE with a white-noise driv- 
ing force of L (t): 

j\(t)+<B ac r\(t)=L(t) (15) 
The autocorrelation function C L is 8-correlated with a strength 

C L (t, t') = <£ (*)£(*')> =*8(f-f') (16) 

Ornstein and Uhlenbeck [H], Doob fl3ll . and van Kampen [0] 
use the integration method to find the statistical properties of 
non-white noise, or colored noise from LE. We, on the other 
hand, solve LE as a I st -order DE. The general solution of a 1^- 
order inhomogeneous DE is a linear superposition of a homo- 
geneous solution (r\h) and a particular solution (r\ p ). Hence, 

^(0=^(0+^(0 

= Ti(0)-exp(-© flc f)+ / ds-exp(-(Oac(t-s))'L(s) 
Jo 

(17) 

From Eqn. (fTTlh the stochastic process at the next time step 
t + At can be obtained. 

V[(t + AO 

= il(0)exp(-© flc (f + At)) 

pt + At 

+ / ds exp(-(d ac (t + At-s))L(s) 
Jo 

/t + At 
ds exp(-G) ac (t + At - s))l(s) 

^ v ' 

!H{t, t + A* ) 

(18) 

Let 9({t,t + Af) be the second term of Eqn. (fT8l) . 



d.s-exp(-(D flC (f + A? - s)) -l(s) 

(19) 

By transforming the variables of integration, we can obtain 

rAt 

#(0, At) = / ds-exp(-Gv(Ar-£)) •£(£ + t) 
Jo 

rAt 

= exp(— Qbac^t) • / J,?- exp(co flC i r ) • L (s + 
Jo 

(20) 

By squaring Eqn. (l20l) , we arrive at 



A. Stochastic Properties 

Langevin Equation governs an O-U process. Hence, if we 
use an O-U process to model GMPS current fluctuations, we 
need to solve LE. By solving the convenient 1^ -order linear 
LE, we were able to extract more statistical properties of the 
O-U process besides Eqns. Q, ([5]), and (©. 

Let us first consider a I st -order SDE of the form of LE. 

m=m)+m an 

Here r\(t) is non- white Gaussian noise with the autocorrela- 
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#- 2 (0, At) 



rAt rAt 

exp(-2to ac At) / / dsdsexp(® ac (s + s))l(s + t)L(s + t) 
Jo Jo 



(21) 



The statistical properties of a random variable can be inves- 
tigated by the calculations of various moments. We calcu- 
late the first and the second central moments 9 by averaging 
Eqns. (f20l) and (I2T1) over an ensemble of particles. The first 
two moments determine the complete statistical properties of 
the O-U noise because it is a zero-mean Gaussian process. For 
zero-mean Gaussian, the 1 st moment vanishes. 



(#(0,Af) ) = 



(22) 



Accordingly, keeping in mind that the O-U process is a sta- 
tionary process, the 2 nd moments boil down to 

(^ 2 (0,A0 ) 

rAt rAt 

= exp(-2a> flC Af) / / dsdsexp(a> ac (s + s))( l(s)l(s] 
Jo Jo 

rAt 

= A exp(— 2(ti ac At) / dsexp(2(ti ac s) 
Jo 



2au 



■{l-exp(-2a> flC Af)} 



(23) 



The second moments of H can be expanded in a closed form 
as in Eqn. (l24l) . 



(#- 2 (0, AO 



2®ac 

■ &At 



1 - exp(-2co< 



c At)} 



r 8 

■[2(© flC Af) - 2(co flC A0 2 + ^(co ac A0 3 - 



1 



2 2 1 3 

%t + 3 V - 3V + 



(24) 

with being co ac Ar. What determines the ( rt 1 ) is G) flC Af, 
which is time step At in units of autocorrelation time x ac , not 
autocorrelation time, or time step by itself. This module is 
designed to generate O-U stochastic noise r\(t) that is to be 
applied to macroparticles in the form of magnetic-field pertur- 
bation: autocorrelation time {l ac ), time step (At), and noise 
strength (A). 



9 When a mean value of a variable is included in the moment calculation, it 
is referred to as central moment. 



B. Box-Muller-Like Transformation 

The Box-Muller (BM) transformation HI, Q1J is intrin- 
sically for generating independent Gaussian white noise — 
which is a limiting case of physical noise — from inde- 
pendent uniform random deviates. In order to generate 
exponentially-driven Gaussian stochastic noise, an exponen- 
tial factor, exp(— coA^) is first multiplied by the stochastic 
noise r\(t) at present time t. Then, a root-mean- square (rms) 
value of 9i (0, At) is added to compute the noise at the next 
time step t + At. 



r\(t + At) =exp(-coA?) 'i\(t)+Cw \J (t, t + At) 2 ) 



■■ exp(-coAf) -r\(t) + Cw y (& (0, A^) 2 ), 



(25) 



where Cw denotes random deviates from a rectangular dis- 
tribution (or white noise). What Eqn. (l25l) implies is that to 
generate r\(t + At), one needs to know T| (t) only. This takes 
advantage of the powerful Markov property of the O-U pro- 
cess in numerical calculations. Upon providing with stochas- 
tic parameters, the variant of the BM transformation is capa- 
ble of generating a wide spectrum of stochastic noise: colored 
noise, non-white noise, off-white noise, etc. Sample paths of 
different noises that are generated from the new noise mod- 
ule are plotted in FIG. HJ the autocorrelation time (z ac ) ranges 
from 10" 3 x7b (7b denotes one revolution period.) to 10 4 x 7b. 
The time step is fixed at one revolution period at the Booster 
injection energy. FIG. [H demonstrates that the autocorrelation 
time governs the pattern of sample path. It is therefore evident 
that the pattern of all sample paths are aperiodic. More details 
of the non- white noise algorithm can be found elsewhere fl7ll . 



C. Application of Noise to Macroparticles 

Since current fluctuations are directly proportional to 
magnetic-field fluctuations, in the noise model the ripple- 
current measurements are translated into magnetic-field fluc- 
tuations as in Eqn. i 



K 



imag 



K 



imag 



AK 



imag 



K 



imag 



■( 



i + 



AK, 



imag 



/K (> 



imag ^ 

(26) 

where imag denotes magnet index for differentiating between 
each individual main magnets. In order to distinguish field 
fluctuations at each type of magnet (F, or 7)), Ki mag is factored 
out, and the amount of field variation (AKi mag ) is normalized 
by Kjmag as a perturbation term. 
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FIG. 4: Sample paths of the Ornstein-Uhlenbeck noise over 1,000 tracking turns; the autocorrelation time (x ac ) ranges from 
10 -3 x To to 10 4 x To, where To denotes one revolution period; the horizontal axis is turn number and the vertical axis 

is noise amplitude. 
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It should be noted that according to experimental measure- 
ments (see section [Vllb . the amount of ripple current (A/) is 
positive above the baseline of a sinusoidal current waveform. 



Hence, absolute values of AKi mag ( 



AK 



imag 



) are taken to rep- 



resent measured ripple current as in Eqn. ([26b . 



VI. GRADIENT-MAGNET POWER-SUPPLY (GMPS) 
SYSTEM 

The Gradient Magnet Power Supply (GMPS) System for 
the Booster synchrotron powers a total of 96 main gradient 
magnets 10 . A resonance system is selected in order to re- 
duce the size and the cost of the power-supply system. The 
Booster magnet system consists of 48 LC-resonant magnet 
cells. A focusing magnet (F), a defocusing magnet (D), a 
choke, and a capacitor bank constitute an individual magnet 
cell. In turn, 48 focusing and 48 defocusing magnets are con- 
nected in series by common buses. Since the gradient mag- 
nets are powered by four independent power supplies (GMPS) 
that are symmetrically inserted in the LC-resonant system, the 
GMPS voltages to ground (V+g and V_g) can be kept as low 
as possible. Each GMPS drives current at the fundamental 
frequency of 15 Hz through a string of 12 magnet cells. The 
GMPS system includes dual three-phase Silicon Controlled 
Rectifier (SCR) bridges connected in series, and fed by a 12- 
phase 13.8-kV bus with shunt (or stray) capacitors connected 
to ground. The components of the Booster GMPS system are 
summarized in Table U 



TABLE I: The Booster Gradient-Magnet System 



Component 


No. 


LC-resonant magnet cells 


48 


focusing magnets 


48 


defocusing magnets 


48 


chokes 


48 


capacitor banks 


48 


GMPS 


4 


gradient magnets / cell 


2 


choke / cell 


1 


magnet cells / GMPS 


12 



VII. NOISE-MEASUREMENT TECHNIQUES AND 
ANALYSIS 

A. Common-Mode Noise and Differential-Mode Noise 

Starting from summer 2005 through winter 2006, we re- 
peatedly conducted direct measurements of 15 -Hz current on 



the main bus line and common-mode and differential-mode 
voltages at each of four GMPS units. From a series of mea- 
surements, we confirmed that ripple current and common- 
mode voltages have consistently been detected and they are 
not of seasonal behavior at each individual GMPS unit. A 
15 -Hz current waveform over 3 cycles is shown in FIG. \5\ 
Frequency and period of the current are shown in the shaded 
boxes on the figure. Two of four GMPS units located in 
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FIG. 5: 15 -Hz current waveform taken at GMPS #1 



the East Booster gallery are pictured in FIG. As illus- 




FIG. 6: Two GMPS Units 

trated by FIG. [71 the waveforms of both V + g and V-g sig- 
nals are sampled at the two leads on the GMPS control rack. 
Utilizing a digital oscilloscope, 11 common-mode voltages and 
differential-mode voltages are calculated as follows: 



V CM =V+G + V_ G 
Vdm = V+G - V-G 



(27) 



The waveforms of V+g and the inverted V-g (V-g) are over- 
laid for easy comparison on the same scale in FIG. [8] In 
addition to FIG. [SJ overlaid CM voltages are plotted against 
V+g and V-g for each GMPS in FIG. Peak-to-peak mea- 
surements quantify the potential differences and cursor-key- 
function determines phase lags between two signals. Potential 



10 The gradient magnet is referred to as the combined-function magnet of 
the Booster. 



11 The model name of the digital oscilloscope used for the measurements is 
Agilent 54622A, and the part number 54622-97014. 
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♦ I cm: Common-Mode Current 

♦ I dm: Differential-Mode Current 




-adient Magnet 
Power Supply 



Shunt Capacitors 



FIG. 7: Common-mode current and differential-mode 
current at the Booster GMPS 



TABLE II: Difference in voltage amplitudes and phase lags 

at each GMPS 



GMPS No. V H 


-G (V) V- 


-g(V) |^| 


AX(ms) 


GMPS 1 


1.577 


1.905 20.8 % 


0.6 


GMPS 2 


3.232 


1.699 47.4% 


4.0 


GMPS 3 


1.598 


1.740 8.9% 


1.4 


GMPS 4 


1.581 


1.743 10.2% 


4.6 



transductor electronics and a current amplifier (TA22 Texas 
Instruments) are used for signal amplification. FFT impulses 
of current waveform of one cycle is displayed in FIG.Q21 The 
vertical scale is 20 dB/div and the horizontal span is 1,670 
Hz. 



differences are displayed in FIG. [lOl In principle, the wave- 
forms of V+G and V-g are supposed to be 180 degrees out of 
phase. However, as displayed in FIG.[TTJ substantial amounts 
of phase lag are found at each GMPS. We calculated the frac- 
tional difference in amplitudes (|AV/V|) by taking the dif- 
ference between V + g and V-g and normalizing the difference 
by V+g- It was discovered that substantial amounts of po- 
tential differences in V + g and V-g are present at each GMPS 
unit. The counterparts of DM voltage are shown in FIG. [121 
The DC output of the power supply is filtered with a L-C net- 
work, with the aid of a 15-Hz low-pass filter, to smooth the 
differential-mode sawtooth waveform at all GMPS units. The 
waveforms of CM voltage (Vcm) are noticeably fast fluctuat- 
ing, which will induce additional current fluctuations in the 
system. Peak-to-peak amplitudes of V+g and V-g are plot- 
ted in FIG. [lOl Phase lags between waveforms of V+g and 
V-g are plotted in FIG. [TT] The voltage divisions are set 
to 500 mV /div and the sweep speed is set to 20 mV / div in 
FIG. [8] through FIG. [TO) From these measurements, it was 
found that the mismatch of potential differences and phase 
lags between V+g and V-g for each GMPS unit are different. 

It was found that the following are the two root causes of 
common-mode noise arising from each GMPS unit: 



(1) additional phase lags (AX) between V+g and V-g 

(2) amplitude (potential) difference between V + g and V-g 



When V+g and V-g are added in a point-by-point fashion on 
the scope, they do not cancel out each other. Instead, the rip- 
ples on each waveform add up and the common-mode voltage 
thus stands out. 

The potential differences and the phase lags measured for 
individual GMPS units are summarized in Table HH Of four 
GMPS units, fractional potential difference in GMPS #2 is 
the largest and the worst. FIG. [13] shows the ripple current 
on a linear ramp of the sinusoidal waveform. The currents 
were sampled directly from the magnet bus line. Since the rip- 
ple currents are such a small fraction of the reference current, 



B. Power Spectral Density of Noise: FFT Analysis 

The measured common-mode voltages from all of four 
GMPS units and the current signal with ripples are Fourier- 
analyzed. To provide real-time proof of the presence of of- 
fending interference in the power-supply system, all the sig- 
nals are analyzed on the fly without being transported to any 
commercial software for the post-measurement analysis. We 
performed real-time analysis with the aid of the built-in FFT- 
function feature on the scope. As FIG. \5\ through FIG. [131 
the real-time graphics were saved on the scope at the time of 
measurements. The resolution of a resonant peak, or FFT bin 
size, is determined by the FFT sampling rate and the num- 
ber of points. The number of points on the scope is fixed 
at 2048, such that the FFT sampling rates and the span of 
the frequency domain are controlled in accordance with the 
Nyquist sampling theorem. In addition, in order to enhance 
spectrum resolution around the frequency peak, the Harming 
window was selected over flat- top, rectangular, and Hamming 
windows. The selected Hanning- window function is given in 
Eqn. d28]): 



1 



cos 



27W 



N 



i) 



(28) 



in which 74^/(7), t and N denote Hanning- window function, 
time, and the number of samples, respectively. 

The following is a list of the FFT settings used on the scope: 

(1) FFT sampling rate, f s = 3.53 kSa/ s 

(2) FFT bin size, Af = 1.04 Hz 

(3) Frequency-domain span = 1.67 kHz 

(4) Horizontal scale = 167 Hz/ div 

(5) Vertical scale = 20 dB/div 

According to the Nyquist sampling theorem, the oscillo- 
scope determines FFT sampling rate from the chosen span of 
frequency domain. 
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GMPS #1 through GMPS #4. As indicated on the 
upper edge of each display, the voltage division is 
set to 500 mV/ div and sweep speed to 20 ms/ div. 
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FIG. 9: The waveforms of Vcm are plotted against those of 
V+G and inverted V_g- Starting from top to bottom, 
each display corresponds to the GMPS #1 through 
the GMPS #4. 
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FIG. 10: Peak-to-peak amplitudes (Pk-Pk(l) and Pk-Pk(2)) and frequencies of V + g and inverted V_g 



11 







> 














JU\j 








\\ 
































































AX = BOOus 

1 /AY - 1 KKKTM 










r 






AY( 1 } - 


3.31V 























Print to: 1 Format ] Factors 
J . TIF image J . JL 



(a)[GMPS #1]: phase lag (AX) is 0.6 ms 




(b)[GMPS #2] phase lag (AX) is 4.0 ms 



















r 




















































a: 


< = 1.40ms 
















1/ 

A 1 


AX = 714.29Hz I 
r- c i i = i _^qv 








; 
























=F 


mnt to aisK me: rK» 
* Print to: I ^ Format 
. Disk J . TIF image 


| Factors J 

L -i 1 r 


i -i i ^ J 



(c)[GMPS #3] phase lag (AX) 1.40 ms 

























- 


























































N. 






















AX = 4.60ms 












- — i - 






itcij = 


1.73V 






















mnt to oisk in 
* Print to^^H 
, Disk J^Tl 


ormat T Factors T T T — 1 
F image J .. _| J . 1" 1 _| I ♦ J 



(d)[GMPS #4] phase lag (AX) 4.60 ms 

FIG. 11: phase lags between V + g and inverted V_g 
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C. Parameterization of GMPS Noise 

For stochastic noise models, the autocorrelation time 
x ac can be viewed as a memory span, or a measure of the 
dependence of the same stochastic values at two distinct 
times it and t'). In this subsection, the measured current 
fluctuations are parameterized with the three stochastic 
parameters introduced in preceding sections HV1 and IVl 

(1) time step (At): The entire Booster magnet system is di- 
vided into four quadrants. Each quadrant made up of a 
string of 24 magnets in series connection is driven by one 
GMPS. Current fluctuations (AI/I) from each GMPS are 
transmitted to all magnets in each quadrant of the ring. As 
such, all of the 24 magnets experience the same amount 
of ripple current at an interval of the time step. Hence, 
the time step, or noise- sampling rate is chosen to be one 
revolution period (7b = 2.2 /us) at injection energy of 400 
MeV. 

(2) autocorrelation time, or correlation time (T flC ): On the 

basis of direct current measurements from a main bus line, 
the ripple currents are repeated above the base current, 
or reference current at an interval of 1.5 ~ 1.7 (ms) (see 
FIG. [13]). Therefore, about the duration of 1.5 ~ 1.7 (ms) 
is chosen to be a proper autocorrelation time for additional 
current fluctuations originated from each GMPS. 

(3) noise strength Based upon the amplitudes of ripple 
current (A// 1) on a linear ramping portion of a sinusoidal 
current waveform (cf. FIG.fTSK the rms value of fractional 

current fluctuation AI/I is on the order of 10 -4 . For 

rms 

verification purpose, histograms of the O-U noise gener- 
ated from each of four noise nodes that are symmetrically 
inserted around the Booster ring are plotted. As shown 
in FIG. [141 the rms values of histograms are on the same 
order as those of measured noise strengths. 
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FIG. 13: ripple current on a linear ramp of the sinusoidal 
waveform. The current are sampled directly from 
the magnet bus line. 
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FIG. 14: Histogram of the amplitudes of noise generated at 
each random noise node. 



The autocorrelation function of a signal, or the power spectra 
can be measured by means of FFT. 

S^) = ^J\- i ^c^x)dx (29) 



e' m S^(03)d03 



(30) 



According to the Wiener- Khinchine theorem [4], spectral 
density is the Fourier transform of the autocorrelation func- 
tion C^(t) = ( ^(x)^(f + x) ) for stationary processes: 

Sp(o) = 2 ( £(t)5(* + x) )cos(cox) dx, (31) 
Jo 

with S^(co) being spectral density of a stochastic process 
In FIG. [151 FFT is performed with one-cycle range of time 
data from 15 -Hz current. For FIG Q3] (a), the horizontal 
scale is 167 Hz/ div, and the vertical scale is 20 dB/div. For 
FIG.[l5](b), the vertical scale is 10 -1 to indicate power attenu- 
ation from 1.0. The power spectral density of the O-U noise is 
closely matched to that of the measured ripple current shown 
in FIG.H3 



D. Equivalent- Circuit Model 

To find out whether there are any offending resonances 
floating around the Booster magnet system, acting as noise 
amplifiers, the equivalent circuit modeling was also carried 
out. The equivalent circuit of one single LC-resonant cell is 
drawn in Figure [T6l 181. Since a string of 24 magnets in a 
quadrant of the Booster magnet system are connected in se- 
ries, they are treated as one transmission line. We employed 
the B 2 SPICE El A/D Version 4, which is one of many ver- 
sions of commercial SPICE simulators. The results of AC 
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FIG. 16: Equivalent circuits of a focusing magnet, a 

defocusing magnet, and a choke that comprise each 
magnet cell 



analysis of the equivalent circuits are shown in FIG. [T7J The 
current is peaked at 15 Hz and a cluster of minor peaks are 
found in a few kHz range. It is speculated that the offend- 
ing resonances above 15 Hz in higher frequency region could 
amplify the power supply noise, when the noise frequencies 
coincide with those of resonances. The presence of the reso- 
nances could augment the formation of beam halo, eventually 
resulting in beam loss during the injection cycle. 



VIII. TRACKING AND BEAM DIAGNOSIS 

A. Simulation Parameters 

A comprehensive set of machine parameters for the Booster 
ring at injection energy is given in Table [TTT1 Some param- 
eters listed in Table (TTTJ are derived from lattice parameters 
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FIG. 17: [SPICE simulation]: Current vs. Frequency: 

current flowing through a string of 12 magnet cells 
driven by one GMPS. Progressing from top to 
bottom, the lines correspond to magnet cell [1] 
through magnet cell [12]. 



specified in the Booster design lattice (version 1.1). Salient 
ORBIT-FNAL simulation parameters including space-charge 
calculations are listed in Table [IVl A round beam with ax- 
isymmetry is first injected into the Booster ring before track- 
ing. This ensures that we can solely investigate the noise ef- 
fects under space charge alone. Optics functions (oc(z), 
y(z), 8(z)) are computed with the Booster design lattice using 
MAD (version 8.23) prior to particle tracking. According to 
the latest measurements and actual machine operation param- 
eters, a careful choice of the other simulation parameters are 
made. 



B. Parallelized Space-Charge Calculations 

A total of 330,000 macroparticles were tracked for the full 
injection cycle of about 2.2 (ms). A grid of 64 x 64 cells was 
used for transverse space-charge calculations, and 32 bins for 
longitudinal space-charge calculation in the Particle-In-Cell 
(PIC) space charge model. During the course of tracking, a 
total of 809 space-charge kicks were applied per revolution. 
This corresponds to about 2 kicks /m, or 58.6 cm /kick. In 
terms of betatron oscillations, about 121 kicks per horizon- 
tal betatron oscillations, and about 119 kicks per vertical be- 
tatron oscillations, and about 17 kicks per magnet cell. In- 
cluding both longitudinal and transverse space-charge calcu- 
lations, each parallelized calculation required about 8 hours 
on forty-eight 2.0-GHz worker nodes. With space-charge bin 
numbers of (64 x 64) x 32 fixed, rms emittances from track- 
ing different numbers of macroparticles are calculated. As 
illustrated by FIG. Q13 when the total number of macroparti- 
cles amounts to above 330,000 after injection is complete, the 
time evolution of rms emittance converges with stability. As 



TABLE III: Machine Parameters of Fermilab's Booster at 

Injection Energy 



such, considering practical computing time and the number of 
macroparticles assigned to each space-charge bin, we deter- 
mine that the number of 330,000 macroparticles is sufficiently 
large for accuracy. Each macroparticle in the noise model rep- 
resents O(10 5 ) real particles, or protons in the Booster. Addi- 
tional beam diagnostic calculations, such as invariant action 



ring radius ((R)) 


75.47 (m) 


ring circumference 


474.2 (m) 


injection kinetic energy 


400 (MeV) 


injection momentum 


954.263 (MeV/c) 


synchronous energy (E s ) 


1.328 (GeV) 


(3 (Lorentz factor) 


0.7131 


y (Lorentz factor) 


1.426 


revolution period (7b) 


2.2 (jus) 


revolution frequency (Jo) 


454.5 (kHz) 


no. of injection turns 


11 


injection period 


24.2 (jus) 


cycle time 


66.7 (ms) 


y tr (transition gamma) 


5.4696 


(Xi (momentum compaction factor) 


0.0172 


phase- slip factor ( \r\ ) 


0.458 


&tr, 95, n (95 %, normalized) 


12.0 (7i-mm-mrad) 


RF range 


Jo.lo ~ jz.oJ (MHz) 


Vxo/VyO (bare tunes) 


6.7 / 6.8 


betatron frequency (f^, fc y ) 


318.2/ 363.6 (kHz) 


Q s (synchrotron tune) 


1.147 xl0~ 3 


Cl s (synchrotron frequency) 


3.28 (kHz) 


T s (synchrotron period) 


305 (jjs) 


G z (rms bunch length) 


l.U (mj 


\\ ( \C\X\ CTltllHlTI C\\ r\(*\c\ IllTir'tinTl 1 
jj£ ^IVJllglLUU-llldl UCld lUlldlUll^ 


^ v 1 f) 4 (m^ 


o/ ^lUllgllUU.llld.1 CllllLLallUCy 


w.z, j \p y oj 


DdLcii liiieiisiiy 


^ ozl v 1 n 12 

J.L/H- A 1U 


average beam current (at injection) 


A Ofl A \ 


effective beam radius 


U.UJZj (JTi) 


effective beam-pipe radius 


U.UOJJ (JOCi) 


bunching factor (Bf) 




Z_W ^LUIie SllllLSy 




AP 1 

Po \max 


± 0.15 % 


3.0 xl0~ 4 




33.7 / 20.5 (m) 


Dx,max 1 Dy max 


3.2 / 0.0 (m) 


cell type 


FOFDOOD 


cell length 


20.62 (m) 


gradient magnets / cell 


4 


total gradient magnets 


96 


V r f, inj (RF voltage at injection) 


205.0 (kV/Turn) 


phase advance / cell 


96 (deg) 


po (defocusing bending radius) 


48.034100 (m) 


pF (focusing bending radius) 


40.847086 (m) 
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TABLE IV: Salient Booster Simulation Parameters with 

ORBIT-FNAL 



no. of injection turns 


11 


no. of maximum macrop articles 


330,000 


harmonic no. 


84 


beam kinetic energy 


400.0 (MeV) 


beam intensity (per RF bucket) 


6.0 x 10 10 


transverse beam distribution 


bi-Gaussian 


ring circumference 


474.2 (m) 


Px, in j | in j 


6.274/ 19.312 (m) 


in i C&y in / 


-0.122/0.024 


D x , o Dy t o 


2.581 / 0.0 (m) 


*0, in j yo, in j 


0.0 / 0.0 (mm) 


Eoffset 


0.0 (GeV) 


1 ^kinetic 


5.1 x 10~ 4 


&x, rms, inj \ &y, rms, inj 


1.76/ 1.76 (7l-mm-mrad) 


V r f (RF voltage) 


205.0 (kV/Turn) 


Rwall/Rbeam (for geometric factor) 


2.0 


longitudinal SC bin no. 


32 


transverse SC bin no. 


64x64 


smoothing parameter 


- 10" 6 


no. of total tracking turns 


1,000 



a LSC stands for longitudinal space charge. 
^TSC stands for transverse space charge. 



calculations, were implemented in parallel mode. 
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FIG. 18: Calculations of rms emittances with a varying 
number of total macroparticles and fixed 
space-charge bin numbers; progressing from top to 
bottom, each trace corresponds with 1 1 ,000, 
33,000, 110,000, 330,000, and 1,100,000 
macroparticles in total. 



C. Moments 

In the following subsections, we will present how beam di- 
agnostic quantities are defined and computed for the stochas- 
tic noise model, we define moments which characterize prob- 
ability distributions of a beam, or macroparticles. Since it is 
necessary to consider beam centroids ((x) and (y)) in calcula- 
tions, ORBIT-FNAL employs central moments: 



Ax r = x r — (x r ) 
Ay r =y r - (y r ) 



(32) 



where x r (z) and y r (z) denote real-space coordinates. Be- 
cause of vanishing central moment calculation, beam cen- 
troids themselves ((x) and ( y )) are used for the 1 st moment 
calculations. It is assumed that the density profiles of an actual 
beam in transverse planes are bi-Gaussian. We first injected 
a herd of macroparticles of bi-Gaussian distribution. Then, 
rms beam sizes (c x ^ oiva ^i^nW^ri tVi^ ond 
calculation: 



o y ) are calculated from the 2 nd moment 



1 moments 



2 nd moments < 



X r ) 
{(Vr) 



= ( (Ax r 
= ( (A>v 



> 2 > 
• 2 > 



(33) 



(34) 



The rms beam sizes are important for space-charge study. 
Starting with (bi-) Gaussian charge distribution p(r), we can 
derive transverse space-charge force using Gauss' law and 
Ampere's law. As given in Eqn. d35ll . the transverse rms 
beam sizes (07) determine the range of linear transverse space- 
charge forces (F sc {r)y. 



?sc(r) 
r 



Ne 

Ne 1 
2m y 2 rL b 



(35) 



where N, e, £q, a r , y, r, and are the number of particles per 
length, unit charge, permittivity of the vacuum, rms beam size, 
Lorentz factor, radial distance, and bunch length, respectively. 
The transverse space-charge forces grow linearly with trans- 
verse displacements (x, or y), and scale off with displacements 
larger than 2o r . As the evolution of 2 nd moment calculations 
show in presence of full space charge and power-supply noise, 
the rms beam sizes grow steadily. To understand the time evo- 
lution of rms beam sizes, or rms beam widths in transverse 
planes, the 2 nd moments in real physical space are computed. 

The injection transverse coordinates employed in the 
ORBIT-FNAL are defined in physical space as a function of 
azimuthal coordinate z. Hence, the horizontal coordinates in- 
clude the effects of horizontal dispersion (£> x o(z)). On the 
other hand, no dispersion effect is included in the vertical co- 
ordinates because vertical dispersion ( ( D y o(z)) is set to zero in 
accordance with the Booster design lattice. Consequently, the 
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following relations are implicitly reflected in the macroparti- 
cle coordinates and the calculations of transverse rms emit- 
tances: 




(36) 



Transverse rms emittances are defined in (xp, jPr/Po) and 
(yp, &y/Po) phase spaces, following the MAD [12011 conven- 
tion: 



((Axp) 2 ) <A*pA^> 
(A^Ajcp) ((A^ x ) 2 ) 



(42) 



In Eqn. (|36K xp (z) and yp (z) denote betatron coordinates, and 
D x o(z) and Po denote injection horizontal dispersion and de- 
sign momentum, respectively. In a similar fashion, divergence 
angles are computed: 



\y'r{z) =y p (z), 

where <D f x (z) denotes the slope of horizontal dispersion. 



(37) 



D. RMS Emittances 

As the ORBIT-FNAL employs the 2 nd -order central mo- 
ments in the rms emittance calculations, we need to define 
additional quantities below to define rms emittances. 



Av x =<p x -($ x ) 

A?y =$ -($ ) 



(38) 



Once we define the 2 nd -order central moments of each coor- 
dinate in the 6-dimensional space, we define column matrices 

M 2 , x and M 2 , r 



Axp 
Ax' 



M 2 , * 



Ayp' 



(39) 



With the column matrices M 2 above, we can define 2 x 2L- 
matrices in subspaces of trace space: (xp, xp and (yp, yp. In 
each of Z-matrices, the off-diagonal elements are associated 
with the correlation between position and angle, or energy and 
rf phase. 



E(xp, X ') = {M 2 , X .Ml, 



<(A*p) 2 > (Ax^Ax') 
(Ax'pAxp) ((Ax'/) 



(40) 



in which M T denotes a transpose matrix of M . Using the 
2x2 Z-matrix, an unnormalized rms emittance can be defined 
as, 



£x, rms = A/det Z(xp, Xp) 



((Ax p )2)((Ax' p )2)_(AxpA4^ 

correlation term 



(41) 



- A /det Z(xp, fpjc) 



(43) 



in which the transverse momenta (2> x and <P y ) are normalized 
by the design momentum (2>o). As stated earlier, if a unnor- 
malized rms emittance is multiplied by the Lorentz factors 
(Py), it transforms into a normalized rms emittance with no 
momentum dependence. The Eqns. d39b through (l43l) apply 
likewise to vertical and longitudinal planes. 



(P Y ) v /((Axp)2).((Axp2)_ ( AxpA4)2 

(Py)%, rms 

/ \2 



(Py)^((Ayp) 2 )((Ayp2)_ (A ypAy p ) 
^-V((Av) 2 )((A^)-(AyA^P 



(44) 



IX. IMPACT OF GMPS CURRENT FLUCTUATIONS 
WITH SPACE CHARGE 

After inserting a total of four random-noise nodes, pro- 
vided with characteristic stochastic noise parameters, into a 
Booster ring, macroparticles representing the Booster beams 
are tracked over 1,000 turns in the presence of full space 
charge 12 . 



12 In the present context, full space charge is referred to as both transverse 
and longitudinal space charge, or 3-D space charge 
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FIG. 19: transverse rms emittance growths starting from the outset of injection through 1,000 tracking turns; 




FIG. 20: Time evolution of rms emittances in comparison with actions; (a) horizontal rms emittance vs. horizontal action; 

(b) vertical rms emittance vs. vertical action 
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FIG. 21: Time evolution of 2 nd moments in transverse planes; 
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FIG. 22: transverse emittance growths; the noise and space-charge effects in red and the space-charge effects alone in blue. 

The beam intensity is 6 x 10 9 ppb, and 5 x 10 11 protons in total. 
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FIG. 23: Fractional exclusion of macroparticles at a given 
average action. The blue indicates at the I st turn 
and the red indicates after 1,000 turns. The vertical 
axis on the left plot is in linear scale, and the right 
is on logarithmic scale. 
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FIG. 24: The distribution of actions (j x and J y ) at the 1^ turn 
and after 1,000 turns. O-U noise and 3-D 
space-charge effects are included. Action 
distribution at the I st turn is in blue, and action 
distribution after 1,000 turns in red. 

As shown in FIG. QjO the time evolution of transverse rms 
emittances with the O-U noise 13 coupled to the full space- 
charge effects (red) and with the space-charge effects alone 
(blue). The beam intensity per bucket is 6.0 x 10 10 ppb, and 
the batch intensity is 5.0 x 10 12 protons. To estimate the emit- 
tance growth rate, the relative emittance growths (^) 14 are 

calculated starting from the last injection turns (the II th turn) 
through 1,000^ turn, prior to beam acceleration; this corre- 
sponds to the first 2 ms out of one cycle over 66.7 ms (15 
Hz). In the horizontal plane the relative emittance growth is 
about 7.5 %, and in the vertical plane the growth is 9.3 %. 
A total of 330,000 macroparticles, or 30,000 macroparticles 
per each injection turn are simulated and tracked. Upon in- 
cluding O-U noise representing the Booster GMPS noise un- 



der space charge, the process of beam degradation develops, 
and a more noticeable halo formation is found. As a cross- 
check with the rms emittance calculations, we also compute 
average actions at each tracking turn including the noise and 
the full space-charge effects. The rms emittances and av- 
erage actions are overlaid in FIG. [20] for clear comparison. 
The calculations of both rms emittances and actions manifest 
in such a good agreement that beam degradation is substan- 
tially enhanced due to synergistic mechanism between GMPS- 
current fluctuations and space-charge effects. Here, we use 
the term synergistic mechanism meaning that the total effects 
of GMPS noise and space charge are larger than the sum of 
individual effects. The time evolution of rms beam sizes in 
both transverse planes with space charge alone and with noise 
and space charge are illustrated in FIG. [21] When the noise 
is included, the 2 nd moments, which are beam size squared, 
grow faster than in the case for space charge alone. If we 
lower the Booster batch intensity by an order of one magni- 
tude (5 x 10 11 ) from the present operational batch intensity 
under the same conditions, the emittance growths induced by 
the GMPS noise and space-charge effects are not distinguish- 
able from those of noise alone in the absence of space charge 
as shown in FIG. [22] It should be noted that the space-charge 
effect is intensity dependent. Thus, if the beam intensity is 
lowered, so is the space-charge effects. This is a clear signa- 
ture that only when the space-charge effects are substantial, so 
does the GMPS noise have a substantial impact on the Booster 
beam. In addition to the primary beam diagnostic calculations 
of the rms emittances and average actions, we looked into the 
transverse couplings as well in the next subsection. 

A. Couplings 

The computations of the 2 nd -order cross moment (xy)for 
each case are presented in FIG. [25] A marginal amount of 
couplings are introduced due to the full space-charge effects 
(FIG. [25] (b)). When the noise is included alone in the ab- 
sence of the space charge, couplings are somewhat noticeable 
(FIG. [25] (c)). When the noise and the full space-charge ef- 
fects are included, the transverse couplings are substantially 
amplified. We therefore conclude that the noise impact on a 
beam is dependent upon the strengths of the space-charge de- 
focusing forces in the Booster. What FIG. l23l illustrates is the 
percentage of macroparticles that reside outside of a given av- 
erage action including the O-U noise and space charge. The 
blue markers indicate the fraction of excluded macroparticles 
at a given emittance at the I st turn and the red markers at the 
end of tracking after 1,000 turns. 



Hereafter, the O-U noise means the stochastic noise modeled on the 
GMPS noise measurements. 
£o denotes initial emittance, and Ae = |e — £o | 




FIG. 25: 



Transverse couplings in configuration space 
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FIG. [24] compares the distributions of transverse single- 
particle actions (J x and J y ) at the outset of injection and at 
the end of 1,000 turns. It is evident that noise-induced beam 
degradation is enhanced as the time elapses. 
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FIG. 26: The distributions of the magnitudes of the 
2^-order cross moment (( xy )) 



matrix as described in Eqns. (1481) through (|49t that follow: 



^4, xy = 



^4, xz = 



^4, y Z = 



A(|) 



(46) 



-lit 



£(*p, 4' ^ ^p) = (MxyM?) 

E(vp, y p ,8E, <|)) = (^^>, 



(47) 



in which M T denote a transpose matrix of M . 



det 



((Ax p ) 2 ) (Ax p AfP x ) (AxpAjp) (AxpA^) 

(A^Axp) ((AiP x ) 2 ) (A^Ajp) (A<P x Av y ] 

(AypAxp) (AjpAfP*) ((Ayp) 1 ) (AjpA^) 

(A^Axp) (A¥yA¥ x ) (A^Ayp) ((AiPj,) 2 ) 



(48) 

In the same fashion, we can compute 4-dimensional emit- 
tances and couplings on x — z and y — z planes as well. 



det 



<(Ax p ) 2 ) (AxpAfP.) 

(A^Axp) ((A^) 2 ) 

(5£Axp) (SEAV X ) 

(A(j)Ax p ) (A^AfP*) 



(Ax p 8£) (Ax p A(j)} 

{ATM) <A<P,A<|>> 

<(S£) 2 > {8EA® 

(A^dE) <(A(H 2 ) 



(49) 



B. Coupling Magnitude 

In an ideal system the normalized rms emittance remains 
constant. However, nonlinear space-charge effect and cou- 
plings can induce degradation in beam quality. The in- 
creases of the normalized rms emittances indicate that non- 
linear space-charge effect and couplings induced by different 
machine imperfections are present in the Booster. One of the 
great advantages of the realistic accelerator simulation with 
macroparticle tracking is that we can isolate an accelerator 
system condition to narrow down a specific cause of emittance 
growth under investigation. Therefore, in order to look into 
the transverse couplings, we additionally implemented in the 
ORBIT-FNAL new parallelized calculations of 4-dimensional 
transverse emittances and coupling magnitudes. From 

the determinant of 2x2 S-matrix of beam distribution as 
given in Eqn. (l42l) . a squared 2-dimensional rms emittance 
on the horizontal plane can be calculated. 



F = det 

e x, rms ^ L 



((Axp) 2 } (AxpArP*) 
(A^Axp) {(A<P X ) 2 ) 



(45) 



Likewise, for the counterpart on the vertical plane. 
For the computation of 4-dimensional rms emittances and 
couplings, we first define 4-component column matrices 
(5W4, ^, 9^4, yz , and 5^4, xz) on two planes as in Eqn. d46*l ). By 
means of generating 4x4 S-matrices with the colum matrices 
as in Eqn. j47t . we can calculate the determinant of each 4x4 



det 



<(Ajp) 2 > (AjpAfP,) (Ay^E) (Ay$A§) 

(A^Ayp) {{AT y ) 2 ) (A<P y SE) (A^A*) 

(5£Ajp) (§EAP y ) ((S£) 2 ) (8£A(j)) 

(A^)Ayp) (AtytoPy) (A^E) ((A$) 2 ) 



(50) 



Hence, coupling magnitudes between x — y, y — z, and x — z 
can be calculated as follows: 



Ae 4 



£ 2 £ 2 



y^Cxy( (Ax p Ajp),(AxpArP > ,}, (AypArP*), (A¥ x AT y ) 



23 terms 



A£ 



yz 



£ 2 £ 2 
^y ^z 



Y,C yz { <Ay p 8E),<Ay p A4>), (8EA<P y ), <Atf y A*) 



Ae 4 

^xz 



e?-e 2 



C K ( (Axp8£),(ArpA(^), ( 8£AfP x ), (A^) 

(51) 

where Cry(...) denotes coupling terms as a function of 
(AxpAyp), (AxpAiPj,), (AypAa^), and (A<P x A<P y ). Hence, 
Aei, includes all possible combinations of couplings not only 
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between horizontal and vertical positions, but also between 
positions and divergence angles in transverse planes. Further- 
more, the 4-dimensional couplings can be extended to hor- 
izontal and longitudinal planes, and to vertical and longitu- 
dinal planes. What is illustrated by FIG. [25] is the 2 nd -order 
cross moment of transverse positions ((Ax r Ay r )) to look into 
the transverse couplings in the same vertical scale. In the ab- 
sence of space charge and GMPS noise, transverse coupling is 
not observed. In the presence of space charge, the magnitude 
of coupling is slightly increased but still marginal. However, 
with the GMPS noise alone in the absence of space charge, 
the coupling is more noticeable and some perturbation appear 
over 1,000 turns. When the GMPS noise is coupled to the full 
space-charge effects, the coupling is substantially amplified. 
In FIG.[26l the turn-by-turn calculations of the cross moments 
are presented in a form of histogram from which we extract 
statistics. The distributions are slightly dispersed as each in- 
stability (either space charge, or GMPS noise) is individually 
included. When the GMPS noise is applied to macroparticles 
in the presence of space charge, the RMS value is larger than 
that of the noise alone by about a factor of two. As derived 




C. Halo Magnitudes 

The computation of maximum extent of macroparticle co- 
ordinates in a beam at each tracking turn is implemented in 
the Noise module. The Eqn. (l52b includes only physical co- 
ordinates (x and y) of a maximum-displaced macroparticle at 
the location of a random noise node (1J] . We refer it to as halo 
magnitude (Rh, 2)'. 



%H,2 



\fx 2 +y 2 



Max 



(52) 



In Eqn. (l53lh halo magnitude in 4 dimension, (Rh, 4), which 
includes horizontal and vertical positions and angles of a 
maximum-displaced particle is given: 



'+(7p;.x') 2 +(v/, 



'IV/) 2 

(53) 

where (3 X and (3^ are optics functions at the location of a noise 
node. Calculations of two types of halo magnitudes (%h, 2 and 
%h, 4) yield consistent results. FIG. [28l illustrates the evolu- 
tion of halo magnitudes in green and smoothed data in blue. 
Due to the large oscillatory behavior of the halo magnitudes, 
the data is smoothed. The smoothed curve in FIG. [29] shows 
us with clarity a growing pattern of a maximum-displaced 
macroparticle from the physical center of a magnet aperture. 




FIG. 27: The time evolution of 4-dimensional coupling, 



in Eqn. (|5TK the coupling between horizontal and longitudi- 
nal planes are continually growing when the GMPS noise and 
space charge impinge on the Booster beam. In FIG. [27j pro- 
gressing from bottom to top, each trace line corresponds with 
each of the following cases: (1) without space charge, nor 
GMPS noise, (2) GMPS noise alone, (3) space charge alone, 
(4) GMPS noise in the presence of space charge. In accor- 
dance with Eqn. (|5TK the vertical axis is in units of (71-mm- 
mrad) 2 . It is evident from the FIG.[27]that transverse coupling 
is synergistically amplified when the GMPS noise is coupled 
to full space-charge effects in comparison with the other cases. 
From the coupling calculations of Az xy and (Ax r Ay r ), we ob- 
tain consistent results; the space charge amplifies the impact 
of GMPS noise on the Booster beam 



S £0S 400 630 80S 1S00 

Turn No. 

FIG. 28: Halo magnitudes (Rmax)'- noise in the presence of 
the space-charge effects; the blue trace in the 
background indicates smoothed curve with spline 
function. 



X. DISCUSSIONS AND CLOSING REMARKS 

The model presented in this paper is the first-ever 
measurement-based stochastic noise model applied to an ex- 
isting low-y accelerator lattice structure through stage ap- 
proach. Utilizing the state-of-the-art parallel computing tech- 
nique for better accuracy, we successfully incorporated and 
tracked a sufficiently large number of macroparticles with 
FFT 3-D space-charge calculations in a practical amount of 
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FIG. 29: Halo magnitudes (R m ax)' noise along with 

space-charge effects (red) vs. space-charge effects 
alone (blue) 



computing time. At first, the new noise module, which can 
generate a wide spectrum of stochastic noise ranging from 
white noise to colored noise, was seemlessly integrated into 
the existing ORBIT-FNAL. We then followed up with discov- 
ering the presence of a substantial amount of offending rip- 
ple current induced by common-mode voltage in the Booster 
power system. However, on the other hand, the differential- 
mode voltage at each individual GMPS is well- smoothed with 
the aid of a 15-Hz low-pass filter installed in each GMPS unit. 
Moreover, the root causes of the presence of common-mode 
voltage at each of four GMPS units were carefully diagnosed. 
As a result of parameterization of the Booster GMPS noise 
from the ripple-current measurements with time step, autocor- 
relation time, and noise strength, FFT power- spectral densi- 
ties between physical noise and modeled Ornstein-Uhlenbeck 
noise are closely matched. 

The foregoing results from particle tracking, with the in- 
clusion of power- supply noise as perturbation and space 
charge as collective instability, make evident that non-white 
noise originating from power supplies under the influence of 
space charge leads synergistically to an enhancement of beam 
degradation phenomena — emittance growth, halo formation, 
and consequential beam loss — at the injection energy of the 
Booster. As mentioned earlier, our investigations evidenced 
that the adverse effects of ripple current are dependent upon 
the strength of space charge. Therefore, as a relevant side, we 
can propose two approaches to coping with the impact that 
ripple current has on charged-particle beams under the influ- 
ence of space charge. The first is to reduce inherent space 
charge forces themselves. Over the past years, the efforts have 
been made to reduce the space-charge effects in the accelera- 
tor system at Fermilab. For instance, in 1993 Fermilab's pro- 
ton linac was upgraded from a beam kinetic energy of 200 
MeV to 400 MeV by adding more klystron tanks in order to 
reduce the space-charge effects in the Booster. Besides, a dual 
RF system with a proper choice of RF parameters, allows us 
to further reduce space-charge effects in high-intensity pro- 
ton machines by means of maneuvering charge distribution 



in longitudinal direction [6]. Accordingly, attendant beam 
degradation phenomena induced by fluctuating current and 
space charge can be suppressed. The second approach is to 
devise instrumental techniques to cancel out common-mode- 
conducted EMI originating from power supplies. In partic- 
ular, as demonstrated by the simulation of the equivalent- 
circuit model of the magnet system serving as an auxiliary 
model, experimental measurements, or detection of harmful 
high-frequency (HF) resonances residing in the magnet sys- 
tem (cf . FIG. [T7]) need to be pursued. Once the presence of a 
cluster of parasitic HF resonances are confirmed, it is required 
that those resonances be damped out to avoid the amplification 
of the adverse influence of power-supply noise on the Booster 
beam. 

Upon including more realistic and non-uniform charge- 
density distribution [17], the effects of space charge increases, 
so does the impact of current fluctuations on the Booster beam 
accordingly. The modeling methodology presented in this pa- 
per is expected to be well applicable to other synchrotrons, 
or storage rings, in which space-charge effects are of con- 
cern. We therefore speculate that power- supply ripple current 
can induce more prominent development of beam degradation 
process in storage rings of space-charge-dominated regime 
over long period of time. 
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